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Abstract 

Motivated by recent experimental work on magnetic properties of Si-MOSFETs, we report a calculation of magneti¬ 
sation and susceptibility of electrons in an inversion layer, taking into account the co-ordinate dependence of electron 
wave function in the direction perpendicular to the plane. It is assumed that the inversion-layer carriers interact via a 
contact repulsive potential, which is treated at a mean-field level, resulting in a self-consistent change of prohle of the 
wave functions. We hnd that the results differ signihcantly from those obtained in the pure 2DEG case (where no pro¬ 
vision is made for a quantum motion in the transverse direction). Specihcally, the critical value of interaction needed 
to attain the ferromagnetic (Stoner) instability is decreased and the Stoner criterion is therefore relaxed. This leads to 
an increased susceptibility and ultimately to a ferromagnetic transition deep in the high-density metallic regime. In 
the opposite limit of low carrier densities, a phenomenological treatment of the in-plane correlation effects suggests 
a ferromagnetic instability above the metal-insulator transition. Results are discussed in the context of the available 
experimental data. 
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1. Introduction 

Recently, there has been much interest in the un¬ 
usual magnetic and magnetotransport properties of 2- 
dimensional (2D) electron gas (2DEG) at low densi¬ 
ties, as exemplihed by the inversion layers in Silicon 
metal-oxide-semiconductor held effect transistors (Si- 
MOSEETs) mSl. or by GaAs quantum wells. Specif¬ 
ically, indirect magnetisation and susceptibility mea¬ 
surements Hi indicate an increase of low-held mag¬ 
netic susceptibility with lowering the 2-dimensional 
carrier density n toward the critical value ric, corre¬ 
sponding to the metal-insulator transition (MIT). This 
might imply the presence of a ferromagnetic transition 
at or near MIT although the situation remains un¬ 
certain due to experimental difficulties. 

On the theory side, the possibility of ferromagnetism 
in a 2DEG at low density was raised in the numerical 


Email address: Denis.Golosov@biu.ac.il (D. I. Golosov) 


investigationsQQ]. Later, a divergence in susceptibility 
with decreasing n toward ric was reported based on an 
advanced renormalisation-group treatment in Ref. 101. 
On the other hand, numerical studies based on the diffu¬ 
sion Monte Carlo technique^ did not hnd any critical 
behaviour of susceptibility in a low density two-valley 
2DEG. 


The effects of hnite thickness of the inversion layer 
on the quasi-2DEG susceptibility were addressed, e.g., 
in Refs. 10,[13.[HI- In these investigations, the appro¬ 
priate diagrammatic summations |0] or the Monte-Carlo 
numerical results 113 are generalised for the quasi- 
2D case by including the appropriate formfactorslll]. 
whereas for interpreting the experimental hndings in 
Ref. ifnl] . orbital effects of the in-plane magnetic held 
(which can be non-negligible for a sufficiently thick 
layer) were invoked. 


We note that within the general framework of 
delocalised-electron magnetism, an important bench¬ 
mark is provided by Stoner-type mean-held theories. 
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The underlying assumption is the presence of a local re¬ 
pulsive interaction between opposite-spin electrons. In 
the case of a quasi-2DEG, where an important role is 
played by long-range Coulomb correlations, such an as¬ 
sumption does not at the first sight appear too realistic. 
Yet we suggest that it might be worthwhile to study it 
in some detail, especially as a relative simplicity of the 
situation allows for a full mean-field treatment includ¬ 
ing the electrostatic effects. We thus aim at general¬ 
ising the classical Ref. Cl to treat the magnetic re¬ 
sponse of a silicon inversion layer, where in addition to 
the long-range Coulomb interaction (which at the mean- 
field level is included in the effective electrostatic poten¬ 
tial), contact (on-site) repulsion between the carriers is 
also taken into account. We note that in the case of sili¬ 
con, recent density functional theory results lfH suggest 
a rather large value of on-site repulsion Uon-site ~ ieV. 
This is somewhat unexpected for a semiconductor with 
covalent bonds, and perhaps provides an additional mo¬ 
tivation for the present study. The validity of our mean- 
field approach will be discussed elsewhere. 


2. Mean-Field Description of an Inversion Layer 


Following Ref. ||13[] . we consider an n-doped Si in¬ 
version layer, with the value of chemical potential fixed 
at the top of the valence band of the bulk silicon. We 
set Fi, to be our zero of energy, and we further assume 
the presence of an acceptor level at £ = E^, - 0, with the 
volume density of acceptors Na- The bottom of conduc¬ 
tion band bends from its bulk value Ec ~ 1.12eV (at¬ 
tained far away from the surface, where the transverse 
co-ordinate z is large and positive) to a variable value 
Ecs < 0 at the surface (z = 0), which is controlled by 
the electrostatic potential tpgate at the metallic gate. 


0(^ “ 0) — ^i^cs Fc) — ^Qgate + ^gate - (1) 

Here, we assume that the potential 0(z) vanishes deep 
within the bulk (more precisely, at z > Zd, where Zd is 
the depletion layer width); —e is the electron charge, 
C is the capacitance per unit area of the oxide layer, 
and Qgate is the (positive) surface charge density at the 
gate, which exactly compensates the induced charges in 
the semiconductor: Qgate - + ni + NAZd), with 

(ill) being the density of the 2DEG carriers with spin 
up (down). The constant difference in work functions 
is absorbed into (pgate- The value of 0 at z > 0 is found 
from the (self-consistent) Poisson’s equation as 

(piz) ^-(Ec- Eos) - 


47Te 2ne , 

— —(NAZd + «! + Na + 

L Jo 

where e « 11.9 is the static dielectric constant of bulk 
Si. The wave-functions ipa (with a =t> i) of the trans¬ 
verse motion are the ground-state eigenfunctions of the 
effective mean field one-dimensional Hamiltonian (we 
recall that in the underlying quantum mechanical prob¬ 
lem transverse carrier motion separates from the in¬ 
plane one, the latter being free): 

'Ha^Eo - - e^('^) + Un^alplaiz) “ - ( 3 ) 

Efalf/a = Eoo-t/^a ■ ( 4 ) 


Here, my is the transverse effective mass [for Si-(IOO) 
surface, my » 0.916 in the units of free electron mass 
nie], El is the in-plane magnetic field in units of gfis 
(bare g-factor times Bohr magnetone), and cr^ is the 
Pauli matrix. Orbital effects of the in-p lane field, which 
have been discussed elsewhere Ol Ulldll . are not included 
in our Ei. U is the contact repulsion (or an s-wave scat¬ 
tering) potential, which is related to the on-site repul¬ 
sion Uon-site on the underlying discrete lattice roughly 
as t/ = a^Uon-site with a x 5.43 A the lattice period. 
While for the Si-(IOO) surface two equivalent electron 
valleys should be taken into account, here we assume 
that the Hamiltonian is diagonal in the valley index. In 
writing Eqs. sa, we assumed that only the ground 
states Eoff (for a =t> i) can be populated by the carri¬ 
ers (electric quantum limit IoId FI Then the respective 
carrier densities at zero temperature are given by 

Ua = -2vEQaO(-Eoa) , V = , ( 5 ) 

where the factor 2 accounts for the valley degeneracy, 
6 is the Heavyside function, and ~ Q.l9Qme is the 
effective mass of the in-plane motion. The net carrier 
density and magnetisation (in units of fig) are then 


ny- -H , M 


«T 


(6) 


The self-consistent mean field scheme is completed by 
an equation for the depletion layer width Zd, viz.. 

271 

- Naz:^ ^Eo- Eos -(nt^T + nizQ , (7) 


where Za - z'4>^(z)dz is the average value of z for 

Spin-Up or spin-down carriers. 


^ Generalisation for the high-density region where several 
transverse-motion levels are active is straightforward. 
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In the U = 0, H = 0 case using an aw^ate lll 31115 1 


Hz) 


-z exp(-fe/2) 


( 8 ) 


with an appropriate b yields the values of Eq which 
agree quite well with numerical solution of Eqs. 0]- 
0. We find a similar behaviour in the U > 0, H = 0 
unpolarised case, where the equation for b now reads 


4m|| 



127re^ 


NAZd + 1 = 0. (9) 



The addition of the second term on the 1. h. s. results 
in a slight decrease of the value of b (and hence in an 
increase of = Z[ - 31b) in comparison to the non¬ 
interacting case, as can be expected on general grounds. 


Figure 1: Schematic representation of the effect of spin polarisation 
on the transverse-motion wave-functions aXU > 0. While in the unpo¬ 
larised case spin-up and -down wavefunctions coincide (dashed line), 
this is no longer true at M > 0 (solid lines). Dotted line represents the 
electrostatic potential energy, -^^(z). 


3. Magnetic susceptibility 


Conventional Stoner picture in the uniform three- 
dimensional case involves the energy shifts of the two 
spin subbands under the combined effect of H and 
the Hartree field without changing the wave functions. 
When using the ansatz ([8]) in the case of an inversion 
layer in the electric quantum limit, this approach yields 
the following expression for the magnetic susceptibility; 


Xst = 


dM 

'm 


H^O 


jUbv 


( 10 ) 


We note that x depends on the carrier density n, due 
to the presence of an n-dependent quantity b in the de¬ 
nominator. This is in turn due to the fact that t/ is a 
three-dimensional contact interaction, whereas the two- 
dimensional contact repulsion between the 2DEG carri¬ 
ers is in general obtained by integrating over z as 


U2D = U 


/^OO 

I il/]{z)il/l{z)dz . 
Jo 


( 11 ) 


In the unpolarised case, while also using the ansatz dSll, 
one obtains U 2 D = 3Ubll6. Since b increases with n 
11131], so does the susceptibility, Eq. (HOl (see below, 
Fig.ll. 

It is however obvious that for U > 0 the transverse- 
motion wave functions i//aiz) do change when spin po¬ 
larisation is induced. In order to clarify this, we will first 
assume that for each spin projection, only the ground 
state Eoa of z-axis motion is populated with carriers. 
While this assumption is quantitatively adequate only 
for smaller densities below 3 ■ 10^^ cm“^, it allows an 
easy insight into the generic qualitative behaviour which 
persists also in a complete (and more cumbersome) 


multi-level picture needed at higher densities, which we 
will discuss later. 

In the presence of spin polarisation, the Hartree field 
of spin-up (majority) electrons creates a potential bump 
for the spin-down electrons [see Eq. 0], with the re¬ 
sult that spin-down electrons are pushed towards larger 
z. This new spin-down electron distribution in turn cre¬ 
ates a (smaller) bump for the spin-up electrons at these 
larger values of z, as a result pushing the spin-up elec¬ 
trons closer to the surface. This situation is shown 
schematically in Fig. [1] While for simplicity we show 
the electrostatic potential (p(z) as a single dotted line (as 
if it does not change when polarisation arises), in real¬ 
ity it is affected by the change in wave functions [see 
Eq. (|2]l], and the resultant new (p(z) must be fed back 
self-consistently into Eq. ©. The interaction energy 
decreases with decreasing wavefunction overlap, so that 
sufficiently large U would imply, at the mean field level, 
a tendency towards ferromagnetism which (unlike in the 
usual Stoner picture) involves a degree of spatial sepa¬ 
ration of opposite-spin carriers along the transverse di¬ 
rection. 

In order to probe this picture variationally, we write 
instead of Eq. (I8]i: 


il^a(z) = 



yjl + Jablz^ exp(-fe„z/2). 


( 12 ) 


Here the parameters and ba control the behaviour of 
wave function at small and large values of z respectively, 
whereas Aq. is a normalisation coefficient. In the absence 
of polarisation, the coefficients here do not depend on 
the spin index a. The values of b and j can be found by 
minimising the appropriate thermodynamic potential G. 
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It is somewhat simpler to work this out at a fixed total 
charge n + NAZd and hence fixed Ecs [cf. Eq. ([1])], as 
opposed to fixed (pgate (while the actual situation corre¬ 
sponds to fixing (pgate, the two approaches are equivalent 
for M = 0 or for an infinitesimal M needed when calcu¬ 
lating susceptibility). In the unpolarised case, we have 





_ 0 

c 

1 

II 

0 


tp{z)dz + 7;- + nEc + 
8v 

Un^ 

1 A . . 

/-»oo 


1 

1- 

0 

— J i/( (z)^(z)dz - 

cNa 



0 

1 ^(z)dz - - 

(n + NAZd){.Ec ~ Ecs) 


0 


Here, the first four terms are transverse and in-plane ki¬ 
netic energy of the carriers, the band-gap contribution 
and the interaction energy; the next two terms are the 
electrostatic energies of carriers and acceptors [see Eq. 
©I, and the last term corresponds to the appropriate 
choice of thermodynamic potential for the case of fixed 
Ecs as explained above. The values of Zd and n are found 
self-consistently from Eqs. (|7]l and (|5]-l6]l, respectively, 
where in Eq. (|5ll one has to substitute for Eqo the varia¬ 
tional energy value, 

+ Ec + 

0 2m|| dz^ 

j T /~*oo /~*oo 

+ -^ j i/(‘*(z)dz-^ j i/(^(z)tp(z)dz. (14) 

Eurther details can be found in the Appendix. 

We note that even in an unpolarised case at t/ = 0, the 
wavefunction (IT^ represents an improvement in com¬ 
parison to Eq. ill, as indicated by a non-zero value 
of the coefficient y > 0 found from the energy min¬ 
imisation. With increasing t/, the value of y decreases 
but stays positive throughout physically relevant range 
of parameter values. When a polarisation is induced, 
M + Q, the parameters of the wave function (ITSli de¬ 
pend on the spin index. As long as polarisation is small, 
Mjn 1, we write to leading order 

fio-= (1+/3)f7, ja^y + 6. (15) 

Collecting the leading-order corrections to the self- 
consistent electrostatic energies of carriers and accep¬ 
tors, as well as to the variational energies E^a of the 
transverse motion, enables one to evaluate the suscep¬ 
tibility. The resultant expression is somewhat cumber¬ 
some and is given in the Appendix [see Eq. M l. 

Results are presented in Eigl2] (a), where the dashed 
line corresponds to the variational calculation based on 
Eq. ([T2I) . The computed value of x exceeds that of 


the non-interacting case (Pauli susceptibility, s v in 
our notation), and increases with increasing n. One can 
see that the variational calculation somewhat overesti¬ 
mates the tendency toward ferromagnetism, suggesting 
a transition at 1.9 ■ lO'^ cm“^ , whereas the nu¬ 
merical solution of mean-field equations [including the 
Schrodinger equation (HJi; see solid line in Eig. |2](a)] 
yields the critical value n\ 5.4 ■ 10'^ cm“^. Impor¬ 
tantly, the “variational Stoner” approach, using Eqs. (© 
and ([Tol l (dotted line in Eig. |2] (a) ), strongly underes¬ 
timates X misses the transition. Indeed, the ratio 
Xstly slowly increases with density, reaching the value 
of 2 only at n ~ 5 ■ lO''^ cm“2 (this is for U - 10“^^ 
erg/cm^, as used in Eig. [2]i ; on the other hand, we find 
that for;^f 5 , to diverge at n = ni, one must increase the 
value of t/ to t/ 4.0 ■ 10“^^ erg/cm^, i.e., by a factor 
of 4. 

If one neglects the effect of wavefunction change, de¬ 
picted schematically in Eig. [1] the Stoner criterion of 
ferromagnetism is determined precisely by the diver¬ 
gence of;^f 5 ,, Eq. ([Tol l as 

2vU2d >1 (16) 

i.e., takes on a purely-2D form. We see that treating the 
z-axis motion properly (which includes taking into ac¬ 
count interaction-induced wave functions change in the 
presence of spin polarisation) relaxes the Stoner crite¬ 
rion. 

The inset in Eig. |2](a) shows the numerical results for 
the behaviour of spin polarisation degree 2M/n above 
the transition. It is seen that M increases continuously, 
and the transition appears to be second-order. This is an 
important difference from a conventional Stoner theory 
in two dimensions where (in the case of a perfectly flat 
density of state and at a fixed n, and at the level of the 
present treatment) the total energy at the critical point 
does not change as IMjn is varied from 0 to 1, and a 
fully spin-polarised state is stabilised beyond the transi¬ 
tion, resulting in a magnetisation jump. In reality, how¬ 
ever, the opposite-spin transverse wave functions over¬ 
lap decreases with increasing M (cf. Eig. [T]l, hence the 
effective two-dimensional repulsion, Eq. (HB and the 
Hartree field also decrease, resulting in a smooth transi¬ 
tion. 

As explained above, the results plotted in Eig. |2] 
(a) are only indicative, as in reality at those densities 
where the critical behaviour is suggested more then one 
quantum-mechanical level of the z axis motion is popu¬ 
lated. In Eig. |2] (b), we show the results of full numer¬ 
ical calculation of susceptibility. These are obtained by 
solving Eqs. das, generalised for the case of multiple 
active levels. With increasing density, carriers begin to 
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Figure 2: Magnetic susceptibility ^ as a function of inversion-layer carrier density n. (a) Results assuming single populated level for the z-axis 
motion. Solid, dashed, and dotted line correspond respectively to the numerical solution of the mean-field equations, variational calculation based 
on Eq. {\2\ . and the “Stoner” value calculated using Eqs. m and (H. The inset shows the numerical solution result for the degree of spin 
polaiisation IMjn, which differs from zero above the coiTesponding critical value, «j % 5.4 • 10^^ cm“^. The parameter values ai‘e Na = 10^^ 
cm“^ and U = 10“^^ erg/cm^, the latter con'esponding roughly to Uon-site = 4eV. (b) Numerical results of the full multi-level calculation, showing 
transition at no ~ 2.4 • 10^^ cm“^. 


populate the 2nd, 3rd, 4th, and 5th level at« 3.2 ■ 10'^ 

cmr^, n 2.6 ■ 10'^ cm“^, n x 1.5 ■ lO'"^ cm“^, and 
n 4.0 -10'^ cm“^, respectively. This gives rise to the 
vertical steps in susceptiblity, which are clealy visible 
in Fig. |2](b). Aside from these, the overall behaviour 
is qualitatively similar to the single-level results of Fig. 
|2](a), although the actual mean-field critical density for 
the ferromagnetic transition, no ~ 2.4 ■ 10^'^ cm“^, is 
larger than n i . At n = no, we hnd 4 active levels for each 
spin projection. The increase of critical density as com¬ 
pared to the single-level results of Fig. |2] (a) is due to 
a broader spatial spread of higher-level wavefunctions 
(which reduces both the effective 2D interaction and the 
wavefunction change under the effect of magnetic field). 
Note that in Fig. |2] variational results are presented for 
the single-level case only, as the required calculations 
become more cumbersome in the case of multiple occu¬ 
pied levels. 

3.1. Phenomenological treatment of correlations near 
MIT 

So far, we found that the magnetic susceptibility y 
increases with increasing density n, with a transition 
to ferromagnetism deep inside the high-density metallic 
region. An important underlying assumption was that 
the in-plane carrier motion is free, which is no longer 


adequate for smaller densities close to MIT. Indeed, in 
this region the dominant role is played by the in-plane 
Coulomb correlations, as indicated by the measure of 
the interaction strength, r, = mj^e^ f(eh? ^fm), exceed¬ 
ing unity. Very recently, it has been suggested lfl^ that 
the available data at low densities can be described by 
a phenomenological model of a 2D Fermi gas with a 
renormalised in-plane mass: 

m± = m± -. (17) 

n-nc 

The Pauli-like in-plane magnetic susceptibility would 
then be 13 


Xp 


V, 


mj_ 

iTthf 


(18) 


When n approaches from above, xp increases be¬ 
cause of the effective band narrowing, which does not 
necessarily imply a magnetic instability (in agreement 
also with Ref. 13|). 

Since there appears to be no reason to expect that the 
transverse motion of carriers is affected by the correla¬ 
tions, we suggest that the effects of a short-range U can 
be taken into account by substituting v ^ v in our treat¬ 
ment as described above. It is readily seen that, for ex¬ 
ample, the variational Stoner susceptibility (fTOl i would 
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Figure 3: Magnetic susceptibility x divided by the bare density of 
states V [see Eq. 0] at low densities. In-plane caiiier mass is 
renormalised according to Eq. (n). Solid, dashed, and dotted line 
correspond respectively to the numerical solution of the mean-field 
equations, “Stoner” susceptibility calculated using Eqs. 0,(10) with 
V —> V, and Pauli susceptibility ) 18 > . The parameter values are 
Na = 10 *^ cm“^, U = 10 “^^ erg/cm^, and ric = 6 - 10 *° cm“^. 


diverge at n = ns, corresponding to a zero of the de¬ 
nominator. For the purpose of an estimate we may use a 
non-interacting (U - 0) value of b{n) El, yielding an 
equation 


ns 


lu. 

4 ns - He 




eb? 


11 

NAZd + —ns 


1/3 


= 1.(19) 


This suggests that a ferromagnetic state is stabilised at 
nc < n < ns ■ Further refinement can be obtained by 
using a more adequate variational ansatz such as ([12]), or 
by solving the mean-field equations numerically. This is 
illustrated in Fig. [3] 

Indeed, we see that while;t'p, Eq. (ITSl l. does not show 
any divergence above both the variational Stoner re¬ 
sult and the numerical solution suggest a transition (at 
n, ss LOShc for numerical and at ns ~ 1.07nc for varia¬ 
tional Stoner; this is for the value of U used in Fig. [3]). 
While the two follow each other rather closely in this 
regime, the numerical result for x (™d hence for the 
density value at transition) is larger than the variational 
one, mostly due to a pronounced difference between 
the exact transverse wavefunction and the variational 
one [given by Eq. (|8|)] in this regime. As explained 
above, the transition is expected to be second-order. In¬ 
deed, numerical results imply that M{n) is continuous 


at the transition point, n - n», and rapidly increases 
with decreasing density, reaching full spin polarisation 
at nf > nc- Full spin polarisation is then sustained in the 
region nc < n < nf. 

Alternatively, taking into account the accuracy of ex¬ 
perimental data, it appears possible that the actual MIT 
takes place at the point nf, whereas nc, which enters 
Eq. ([TTT i. plays the role of an asymptotically defined pa¬ 
rameter (a point where the effective mass would have 
diverged). 

4. Discussion 

We constructed a Stoner-type theory of magnetism 
for electrons in an inversion layer, addressing both the 
behaviour of the system in the metallic high-density re¬ 
gion and the correlated low-density regime immediately 
above the metal-insulator transition. While we specifi¬ 
cally aimed at describing Si-(IOO) MOSEETs, our re¬ 
sults are expected to be qualitatively relevant for other 
2D electron systems of finite thickness. These general 
conclusions are: (i) At higher densities, proper treat¬ 
ment requires taking into account the wave function 
change under the applied in-plane magnetic field (see 
Fig- (B- This effect leads to an increased susceptibility 
in the paramagnetic state, and enhances the tendency to¬ 
ward ferromagnetism. As a result, the Stoner criterion is 
somewhat relaxed in comparison to a more conventional 
treatment, which in turn may give rise to a ferromag¬ 
netic transition, (ii) Even though the overall behaviour 
at low densities is dominated by long-range Coulomb 
forces, the magnetic susceptibility can be significantly 
affected by the on-site carrier repulsion, which is of 
course always present. 

We close with a brief summary of our results for a 
silicon MOSEET. In the high-density metallic regime, 
when one neglects the in-plane Coulomb correlations 
(main body of Sec. O, the relaxation of the Stoner crite¬ 
rion translates into a lower value of critical density cor¬ 
responding to a ferromagnetic transition (in comparison 
to the conventional Stoner approach). Still, the transi¬ 
tion to ferromagnetism occurs at rather large n, and it 
is far from clear whether the corresponding values can 
be reached in case of Si-MOSEETs. However, the in¬ 
crease of susceptibility with increasing n should perhaps 
be observable at large yet readily accessible densities, 
n > 10'^ cm“^. 

Describing the magnetic response of the system at 
lower densities close to the metal-insulator transition re¬ 
quires taking into account the strongly-correlated nature 
of the 2DEG in this limit. This was carried out at a phe¬ 
nomenological level in sec. [Q and results suggest a 
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ferromagnetic instability above the metal-insulator tran¬ 
sition. This finding contributes to the on-going debate 
in the literature UQ]- 

The mean-field transition to ferromagnetism is found 
to be continuous in all cases, and a sublinear behaviour 
of M{H) at larger H is anticipated (due to decreasing 
wavefunction overlap). Another important and novel 
feature of our approach is the spatial separation (in 
the transverse direction) which arises between opposite- 
spin carriers in the ferromagnetic case or whenever spin 
polarisation is induced (the average values of the trans¬ 
verse coordinate of spin-up and -down electrons then 
differ, zt < Zj,). This issue will be addressed in more 
detail elsewhere. 


5. Acknowledgements 

The author takes pleasure in thanking R. Berkovits, 
K. A. Kikoin, B. D. Laikhtman, S. V. Kravchenko, I. 
Shlimak, and L. D. Shvartsman for enlightening discus¬ 
sions. This work was supported by the Israeli Absorp¬ 
tion Ministry. 


Appendix A. Details of the variational calculation 
based on Eq. dllb 


As a first step, one has to determine the values of 
the wavefunction parameters b and y in the unpolarised 
case by minimising the thermodynamic potential G, Eq. 
(fTjl) . Evaluating all terms in Eqs. (fTSl) and (fT4]i in the 
usual limit of Za Zd and differentiating, we find 




r-i + 2 . 


2m 

1271 e^NAZd^ 3Tie^n 


+ :^bnU(69f - 210^ h- 205) - 


eb 


5l2be 


(435^^-286^ H-555)=0, (A.l) 


4mii 


-2 - 6a^ + 24-ya^ H-(1 -i- 24ya^) H- - 

2y 2y^ 


3 , rr 4/^n/- irvcN 967Te^NAZd 4 

= - bnUa {69^ - 105) H- a + 

256 be 

32 be 


(A.2) 


Here, the quantities E and Q depend only on y, and are 
expressed in terms of sine and cosine integrals: 


E 


1/1 1 1 1 

—— ci—— sin —— - si—— cos —- , , 

Vr \ Vr Vr Vr Vr/ 

1111 

Q = Cl-cos-h SI-sin-; 

Vr YE w Vr 


(A.3) 

(A.4) 


a - {\ + 127 )“'^^ is the normalisation coefficient in Eq. 
(fTST i at m = 0 , and ^ is defined according to ^ = a^(l + 
20y), so that zj - Zi - 3^/b. Eqs. (IA.lllA.2T i must be 
solved numerically, together with Eqs. (|7|i and (|5}0. 

When an infinitesimal in-plane magnetic field H is 
applied, the spin-up and -down wavefunctions begin to 
differ from each other according to Eq. (fTSl l. This effect, 
along with the Zeeman splitting of the spin subbands, 
contributes to the magnetisation. 


M = 


(694-210^4-205)- 


3bnUv , 
2048 


-/3 + {69^- 105)4x 


3bnUv 

X - d+vH 

128 


3bUv 9 

1- —(694-210^+205) 


. (A.5) 


The thermodynamic potential acquires a correction, 
G —> G + AG, which is quadratic in /?, 6, and H. The 
values of fi and 6 can be found by minimising AG, yield¬ 
ing magnetisation M [with the help of Eq. (IA.51 l1. and 
hence the susceptibility x- A straightforward, if some¬ 
what tedious, calculation yields 


A = 


3bUv 9 
1 - ^y^( 694 - 210f + 205) - K 


1-1 


(A. 6 ) 


Here the first two terms correspond to the Stoner suscep¬ 
tibility, computed using the wavefunctions (fTSll of the 
unpolarised case [instead of ®; cf. Eq. (fTOl il. whereas 
the quantity 


K = 


9vb^nU^ 
2 ' 6 (AB - C 2 ) 

4c 


—(694 - 210f + 205)^ + 4a^Bx 
64 


x(69f- 105)2 - —(694 -210^ + 205)(69^- 105) 


accounts for the field-induced change of the spin-up and 
down wavefunctions, as illustrated by Eig. [1] The coef¬ 
ficients A, B, and C are given by 


A = 


X 


h^a^nb^ 
2m\\ 
6a^ 5 

yZ 3 y 3 

3ne^n^ 

Aeb 




r 


•24a^ 


1 

84 


• 72 fl‘* + 288 aV 


(1305^-429)4- 


2 ^ ■ 97ie^NAZdna^ 
eb 


+ — Ubn\95 - 69^)a^, 

64 

h^a^nb^ 1 3ne^n^ 9 

B = —-(E - - + 24 - ——(915f - 2638^ + 

Imn L l^eb 


+ 1215) + 
-375), 


24ne^NAZdn^ 9 


eb 


-U br? {-199 f+ 510^- 
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c = 
o 

'32 


h^a^nb^ 


2m\\ 


a 


2 , ,2 


+ 12a + 2 + 6a - 24a y 


37Te^n^ 

64 eb 


(435^-751)a*- 


96ne^NAZdna^ 

eb 


+ — Ubn^a^. 
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